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Abstract. Using data taken by SELEX during the 1996-1997 fixed target run at Fermilab, we study the 
production of charmed hadrons on copper and carbon targets with S~ , p, 7r~, and tt^ beams. Parametrizing 
the dependence of the inclusive production cross section on the atomic number A as A" , we determine a for 
, , D , (2010) , A'^ , and their respective anti-particles, as a function of their transverse momentum 
Pt and scaled longitudinal momentum x f . Within our statistics there is no dependence of a on for any 
charm species for the interval 0.1 < xf < 1-0. The average value of a for charm production by pion beams 
is amoson = 0.850 ± 0.028. This is somewhat larger than the corresponding average tibaryon = 0.755 ± 0.016 
for charm production by baryon beams {S~ , p). 



PACS. 13.85.Ni - 14.65.Dw - 24.85.+p 
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1 Introduction 

The inclusive production of various outgoing hadrons from 
a specific incident beam particle interacting in a nuclear 
target has been studied for many years. The usual char- 
acterization of the process is that in a complex nucleus 
of atomic number A, the single nucleon production cross 
section is increased by a factor A" . If the outgoing particle 
is absorbed with the same cross section as the incoming 
beam particle, then a ~ 2/3. Conversely, production of 
heavy quarks by hard scattering with slow hadronization 
may lead to minimal absorption and give a ^ 1. 

Experiments have used a range of techniques to deter- 
mine a for a wide selection of charm and strange hadrons. 
Some experiments have studied generic charm hadron pro- 
duction by detecting only muons from semi-leptonic de- 
cays [UEIE] . Others have measured a for hidden charm in 
J/ tp production [4l[5l[6] . A third group have looked at spe- 
cific charm or strange final state hadrons [7l[8l[9l lT0l[TTlfT2l 
fT3 l fT4 | [T5 p i6 | [T7]. Until this experiment, no single experi- 
ment covered the whole family of charm hadrons over a 
wide range of kinematic variables pt and xp, using tt"*", 
TT~ , p, and X!~ beams. For the first time these new data 
allow one to look for systematic variations of a for differ- 
ent kinematic regions, different groups of outgoing charm 
hadrons, or different beam hadrons. This broad data set 
is important for comparing to the variety of theoretical 
models that have been proposed. 

The production of strange particles with a proton beam 
shows [7l|8] a strong dependence of a on both xp and 
Pt, with a < 0.5 in some bins of {xF,Pt)- Open charm 
production was measured [lllinilllllllilllllllSlIIllIZl for 
pion, proton, and neutron beams, in different ranges oixp, 
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for D*, _D", £>*, and Ds mesons as well as the but 
most experiments only had one beam particle and/or final 
state hadron, or averaged over several mesons, publishing 
just one value for a. In addition to the previously men- 
tioned muon measurements [Tl[11131lllini[3] , prompt neutrino 
production assumed to be from charm decays is reported 
in |18j . As summarized in |19p20j. many measurements 
of open charm production concentrate in the central re- 
gion xp ^ Q and find a ^ \ for the specific final states 
measured. For charmonium production or muon-triggered 
inclusive charm production at large xp, a decreases, ap- 
proaching 2/3. 

Studying the basic production and suppression mech- 
anisms in charm hadroproduction is important for under- 
standing non-perturbative aspects of heavy hadron pro- 
duction. These results also impact other fields like Heavy- 
Ion collisions (see a review [21]) and Cosmic Ray Physics, 
where air shower Monte Carlo simulations have to take 
into account the production of charm particles. 

We present in this letter a new measurement for a in 
the range of 0.1 < xp < 1, for 14 different open charm par- 
ticles and decay modes, produced by four different beam 
particles. In any single beam / charm-particle combina- 
tion this experiment has statistics that are similar to or 
better than those from any other previous measurement. 

2 Experimental Apparatus 

The experimental setup of the SELEX experiment is de- 
scribed elsewhere We point out the most important 
features of the setup used in this analysis. SELEX is a 
3-stage magnetic spectrometer, designed for high accep- 
tance forward [xp > 0.1) interactions. 600GeV/c nega- 
tive (~ 50% ~ 50% TT") and 540GeV/c positive 
beam particles (~ 92 % p, ~ 8 % tt"*"), interact in five tar- 
get foils, described in Tabled] The physical properties of 
the target foils were measured before the installation into 
the experimental setup, and the thicknesses and positions 
were verified by measuring the positions of the primary 
vertices. Every beam particle is individually tagged by a 
Transition Radiation Detector, and the meson (baryon) 
contamination in the baryon (meson) beam is below 1 %. 

The spectrometer had silicon strip detectors to mea- 
sure the incoming beam and outgoing tracks. Momenta 
of particles deflected by the analyzing magnets were mea- 
sured by a system of proportional wire chambers (PWCs), 
drift chambers and silicon strip detectors. Momentum res- 
olution for a typical lOOGeV/c track was <Jp/p « 0.5%. 
Charged particle identification was performed with a Ring 
Imaging Cherenkov detector (RICH) [23j, which distin- 
guished from TT* up to 165 GeV/c. The proton iden- 
tification efficiency was > 95 % above proton threshold 
(w 90GeV/c). For pions reaching the RICH detector, the 
total mis-identification probability due to all sources of 
confusion was < 4%. 

Interactions in the five target foils were selected by 
a scintillator trigger. The trigger for charm required at 
least four charged tracks downstream of the targets as in- 
dicated by an interaction counter (ICl, IC2. see Fig. [1]), 
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Name 



Material 



Thickness L Position A Density p Air 







[cm] 


[cm] 




[g/cm J 


[/OJ 


S4 


Scintillator 


0.158 


7.27 




1.03 


0.20 


6 


Copper 


0.159 


6.13 


63.5 


8.96 


1.06 


7 


Copper 


0.119 


4.62 


63.5 


8.96 


0.76 


8 


Diamond 


0.220 


3.10 


12 


3.25 


0.82 


9 


Diamond 


0.220 


1.61 


12 


3.25 


0.82 


10 


Diamond 


0.220 


0.11 


12 


3.25 


0.82 


ICl 


Scintillator 


0.200 


2.46 




1.03 


0.25 


IC2 


Scintillator 


0.200 


2.97 




1.03 


0.25 



Table 1. Physical Properties of 
Materials in the Charm Production 
Targets region. The layout is shown 
in Fig. [T] 



no signal in a veto counter (S4) upstream of the targets, 
and at least 2 hits in a scintillator hodoscope after the 
second analyzing magnet. It accepted about 1/3 of all in- 
elastic interactions. Triggered events were further tested 
in an on-line computational filter based on downstream 
tracking and particle identification information. The on- 
line filter selected events that had evidence of a secondary 
vertex from tracks completely reconstructed using the for- 
ward PWC spectrometer and the vertex silicon. This filter 
reduced the data size by a factor of nearly 8 at a cost of 
about a factor of 2 in charm yield. From a total of 15.2-10^ 
interactions during the 1996-1997 fixed target run about 
10^ events were written to tape. 



A", we obtain 



IC1 IC2 



-2 



4 6 
z position [cm] 



Fig. 1. Scale drawing of the charm production target region. 
In addition to the five targets (2 Copper, 3 Diamond) we also 
indicate the location of some of the scintillators used in the 
trigger (S4, ICl, IC2) and the first two planes (Ix, ly, from 20 
in total) of the silicon strip detectors. The physical properties 
of the elements are shown in Table [1] 



3 Data Analysis 



To determine the charm production cross section depen- 
dence on the nuclear mass A, we have to determine the 
number of charm particles produced in any single target, 
and take into account the number of nuclei in the Car- 
bon and Copper targets. Parametrizing the cross section 
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In^ 
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In 4^ 

Ac 



In 



PC Lc Acu 
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Lcu, 

Ncn 



with atomic masses Aq, Aq^, the thicknesses Lq, 
and densities pc, Pcu as shown in Table [TJ and Nq, 
being the number of acceptance corrected events observed 
in the different target materials. This expression for a em- 
phasizes the measurement issues. The first term depends 
on the corrected number of events from each target. This 
means that the acceptance and trigger efficiency for each 
different target have to be understood as a function of 
the kinematic variables xp and pt- The second term il- 
lustrates the requirement for precision in establishing the 
parameters of each target used in the measurement. Any 
uncertainties in density or thickness translate directly into 
uncertainty or systematic shift in a. 

In this analysis, we reconstructed completely charm 



particles in specific decay modes. For 
K-TT+n+TT-. D+ K-n+TT+, Df R- 



At 



pK' 



and the corresponding charge- conjugated modes. 



we used cuts similar to those in previous publications [241 
I25p26|. Secondary vertex reconstruction was attempted 
when the per degree of freedom for the fit of the ensem- 
ble of charged tracks to a single primary vertex exceeded 4. 
All combinations of tracks were formed for secondary ver- 
tices (Xsec < 5) and tested against a reconstruction ta- 
ble that specified selection criteria for each charm decay 
mode. Secondary vertices which occurred inside the vol- 
ume of a target were rejected. The resolution of the pri- 
mary vertex position is on average better than 300 /xm (de- 
pending slightly on the target foil), less than the thickness 
of the target foils and much less than the spacing between 
foils. This permits an unambiguous assignment of the in- 
teraction to a specific target foil. Additional identification 
criteria for the different decay modes required that proton 
and kaon candidate tracks were identified by the RICH de- 
tector to be at least as likely as a pion. Additionally, in the 
case of K^K~Tr^ for the kaon tracks the kaon hy- 

pothesis had to be more likely than the proton hypothesis. 
If a pion candidate track reached the RICH detector, we 
applied as a loose requirement that it had to have a likeli- 
hood of at least 10 %; if the track failed to reach the RICH, 
the candidate was called a pion. The separation between 
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Beam Particle 







Decay Mode 




— 


IT 


P 


+ 

TT 


1 




K-'K+ 




1176 ± 38 


411 ± 22 


245 ± 16 


29 ±7 


2 


DO 


K+n- 




1740 ± 52 


452 ± 23 


437 ± 24 


39 ±7 


3 


D° 


— > n^n^n^ 




1282 ± 50 


467 ± 26 


252 ± 18 


47 ±6 


4 


DO 






1650 ± 60 


488 ± 29 


331 ± 26 


73 ±9 


5 


D+ 






1352 ± 46 


361 ± 23 


248 ± 20 


42 ± 7 


D 


D- 






ono/1 -1- 
zUz4 zt oo 


000 zt Z i 


ooo zt ZZ 


OD zt y 


7 


D*+ D°{K-TT^ 


)n+ 


165 ± 13 


48 ±7 


33 ±7 




8 


D* 


- DO{K+TT- 


)n- 


331 ± 20 


70 ±8 


65 ±8 




9 


D*+ D°{K-n^ 


TT^TT )7r^ 


235 ± 15 


61 ±9 


58 ±9 




10 


D* 


- DO{K+TT- 




446 ± 21 


116dz 11 


80 ±10 




11 


Dt 






118 ± 17 


62 dz 11 






12 


D- 


K+R-n- 




379 ± 26 


91 ± 12 






13 


At 






1130 ±39 


172 ± 15 


240 ± 16 




14 


Ac 


— > pK^TT^ 




313 ± 34 


95 ±13 


42 ± 9 





Table 2. Raw yields (before apply- 
ing any acceptance corrections) for 
the charm particles and modes, for 
the different beam particles, used 
in this analysis. These yields were 
obtained fitting a Gaussian and a 
polynomial representing the back- 
ground to the invariant mass dis- 
tributions. 



the primary and secondary vertices had to be greater than 
eight times its uncertainty, and the uncertainty itself less 
than 0.17 cm; the reconstructed charm momentum vector 
had to point back to the primary vertex, and two of the 
daughter tracks had to have a miss distance with respect 
to the primary vertex of more than V6 times its uncer- 
tainty. For D* states decaying into D^n^ , we required a 
reconstructed D° within ±36MeV/c^ (± 3 times the reso- 
lution) of the nominal mass, and an additional pion from 
the primary vertex. The approximate total yields for the 
different modes and beam particles are shown in Table [2l 
The invariant mass distributions were divided into groups 
for the primary interaction happening in one of the tar- 
get foils, and further in different x^-bins, and in some 
cases also in bins of . We used the sideband-subtraction 
technique to remove the background from the mass distri- 
butions. The resultant 736 different yields are the primary 
data for measuring a. 

The total acceptance (geometrical acceptance and re- 
construction efficiencies) for the different decay modes of 
interest was estimated by embedding Monte Carlo charm 
decay tracks into data events. Events were generated with 
independently-parametrized transverse and longitudinal 
momentum distributions, tuned to match the shapes of 
the data. Detector hits, including resolution and multiple 
Coulomb scattering smearing effects, produced by these 
embedded tracks were folded into arrays of hits from real 
events. The new ensemble of hits was passed through the 
SELEX off-line software. We verified that the acceptance 
and reconstruction efficiency does not depend on the mul- 
tiplicity of the underlying event; additionally we found 
that the multiplicity distributions are nearly independent 
of the target material, as shown in Fig. [21 The acceptance 
is the ratio of the number of reconstructed events over the 
number of embedded events in a particular mode for a spe- 
cific target foil and bin in xp and p^ . As seen from equa- 
tion [1] most acceptance, trigger and filter effects cancel 
in this measurement, only differences in the acceptances 
between the target foils are important; the largest effects 
depend on the lifetime of the different states. For example. 

















Cu <14.3> 

C <13.7> 


1 1 1 1 1 p4JT^~i~^~| 


1 




1 1 1 1 1 1 1 


1 





5 10 15 20 25 30 

Multiplicity 



Fig. 2. Multiplicity distributions (number of tracks in the ver- 
tex region) for events with charm candidates used in this anal- 
ysis, for copper (full red) and carbon (dashed blue) targets. 
The corresponding averages are also shown. 

for yl+ decays in the xp interval 0.4-0.6 the acceptance 
varies from 21.4% to 23.5% between target 7 and tar- 
get 10, while for D^ decays the variation is from 41.4% 
to 31.8 %. We verified our acceptance corrections by com- 
paring the corrected event yields as functions of xp and 
Pt for our three identical diamond targets and found good 
agreement within our statistics. This study was performed 
for all the charm decay modes reported here, as well as for 
the high statistics sample yl° — > pn" . 

For the determination of a we use the number of ob- 
served events, corrected for acceptances, from the three 
diamond targets, but only from the second copper tar- 
get (target number 7). Some fraction of the events with 



SELEX Collaboration, A. Blanco-Covarrubias, et al.: Nuclear Dependence of Charm Production 



5 



Beam 


Mode 


Q 


a 


a 


a 






0.1 <XF < 0.2 


0.2 <XF < 0.4 


0.4 <XF < 0.6 


XF > 0.6 


E- 


1 


0.75 ±0.07 


0.72 ± 0.07 


0.48 ± 0.25 






2 


0.80 ± 0.05 


0.70 ± 0.06 


0.98 ± 0.18 


0.71 ± 1.54 




3 


0.52 ± 0.18 


0.66 ± 0.09 


0.57 ± 0.22 


0.67 ± 0.68 


E- 


4 


0.47 ±0.19 


0.67 ± 0.09 


0.80 ± 0.17 


1.23 ± 0.92 


E' 


5 


0.75 ± 0.09 


0.68 ± 0.07 


0.33 ±0.27 




E- 


6 


0.80 ±0.08 


0.79 ± 0.06 


0.74 ±0.13 


0.84 ±0.59 


E- 


7 


0.86 ±0.24 


0.89 ± 0.15 


0.57 ±0.31 




E- 


8 


0.63 ±0.19 


0.73 ± 0.11 


0.74 ±0.20 




E' 


9 


0.43 ± 0.45 


0.41 ± 0.17 


0.88 ±0.17 




E' 


10 


0.80 ±0.21 


0.80 ± 0.10 


0.84 ±0.14 


0.47 ±0.50 


E- 


11 


1.10 ±0.38 


1.07 ±0.19 






E- 


12 


0.99 ±0.35 


0.79 ± 0.12 


0.87 ±0.16 




E- 


13 


0.70 ±0.20 


0.95 ± 0.08 


0.90 ±0.10 


0.83 ±0.17 


E' 


14 


1.32 ±0.25 


0.74 ± 0.24 






TT^ 


1 


0.86 ± 0.14 


0.82 ± 0.11 


0.25 ±0.24 






2 


0.89 ± 0.12 


0.78 ± 0.11 


0.96 ±0.17 




7r~ 


3 


0.85 ±0.32 


0.87 ± 0.12 


0.82 ±0.17 


1.00 ± 0.18 


7r~ 


4 


1.04 ±0.22 


0.88 ± 0.15 


0.59 ± 0.19 


0.75 ± 0.29 


7r~ 


5 


0.37 ±0.34 


0.70 ± 0.13 


0.79 ±0.17 


1.27 ±0.37 


7r~ 


6 


0.76 ±0.18 


0.75 ±0.10 


0.89 ±0.16 


0.83 ±0.21 


7r~ 


7 


1.50 ±0.56 


1.17 ±0.33 


0.68 ±0.36 




7r~ 


8 


0.34 ± 1.06 


1.11 ± 0.24 


0.77 ± 0.32 




7r~ 


9 


1.54 ± 0.60 


0.95 ± 0.25 


0.29 ± 0.46 


1.00 ± 0.38 


7r~ 


10 




0.99 ± 0.22 


0.70 ± 0.25 


1.01 ± 0.22 


7r~ 


11 




0.42 ± 0.49 


0.83 ± 0.58 




7r~ 


12 




0.87 ± 0.24 


0.68 ± 0.44 


1.52 ± 1.04 


TT 


13 


1.42 ± 0.54 


1 08 ± 18 


84 ± 27 


80 ± 37 


TT 


14 






95 ± 37 




p 


I 


KJ'OXJ -1— XJ.A-O 








p 


2 


f) 77 + 1 9 


n C7 _(_ n 1 o 


f) Qo _|_ n 41 






3 


n 77 -1- f) Q4 


n KQ _(_ n Of) 


f) 74 -L n 07 






4 


0.61 ± 0.35 


0.45 ± 0.20 


0.37 ± 0.54 


0.93 ± 2.32 


7} 


5 


0.50 ± 0.26 


0.71 ± 0.14 






7} 


6 


0.94 ± 0.18 


0.80 ± 0.12 


1.02 ± 0.26 




7} 


7 




0.61 ± 0.35 


0.47 ± 0.66 




v 


8 




1.03 ± 0.19 


0.35 ± 0.65 




v 


9 






0.48 ± 0.48 






10 




0.73 ± 0.23 


0.22 ± 0.46 




P 


13 


0.51 ±0.77 


0.44 ±0.23 


0.79 ±0.22 


1.03 ±0.30 


P 


14 




1.08 ±0.80 


0.78 ±0.55 




7r+ 


1 


0.37 ±0.65 


1.08 ±0.29 


0.23 ±0.65 




7T+ 


2 


0.68 ± 0.34 


1.23 ± 0.22 


0.38 ±0.68 




7T+ 


3 




1.00 ± 0.41 


0.91 ±0.35 


1.15 ±0.88 


77+ 


4 




0.74 ±0.32 


1.57 ±0.51 




7r+ 


5 


1.11 ±0.67 


0.43 ± 0.38 


0.58 ±0.80 




1T+ 


6 


0.59 ±0.76 


0.90 ± 0.25 


1.28 ±0.33 


1.03 ±0.60 



Table 3. a- values for the different 
charm particles and decay modes, 
for different beam particles. Only 
statistical uncertainties are shown. 
We do not determine a for a spe- 
cific mode and/or bin of xf if the 
number of observed events (after 
sideband-subtraction) is below 1. 



the interaction in target number 6 where vetoed in the 
hardware trigger by the "S4" scintillation counter due to 
back-splash from the interaction. In the attempt to correct 
for this vetoing we encountered systematic biases which 
would increase the combined statistical and systematic 
uncertainties more than if wc ignore completely the in- 
teractions from the first copper target foil. Nevertheless, 
the results using target 6 are consistent with those without 
it. 



4 Results 

For any single mode, for four different beam particles, we 
calculated a according to equation [TJ The results are pre- 
sented in Table [3l Since any single measurement has a 
large statistical uncertainty, we average our results, weighted 
by the individual uncertainties, for different groups of in- 
terest. 

In Fig. El^a) we show a as a function of xp for all data, 
i.e., averaged over all charm and anti-charm modes and 
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Fig. 3. Average a as function of 2:f for all observed final states (a) and for charm and anti-charm (b). The data points are 
slightly offset to avoid overlapping of the error bars. Reference a values of 2/3 and 1 are shown as dotted lines. The points at 
xf > 1 show the average assuming that a does not depend on xf- 



all beam particles. In Fig. [3{b) we separate the charm 
and anti-charm final states and show a averaged over all 
decay modes and beam particles for each type of charm 
quark. 

In Fig. m^a) we display the dependence of the average a 
on the beam particle type: meson or baryon. All charm and 
anti-charm decay modes are averaged. In Fig. [4l[b,c) we 
separate the charm or anti-charm decays into leading and 
non-leading classes. Recall that leading charm processes 
are those in which the produced charm hadron carries at 
least one valence quark of the beam particle. Non-leading 
charm processes have no valence quarks in common be- 
tween the beam and charm hadrons. In Fig. |4ljc) we ex- 
clude the and C as these states could be produced 
preferably via excited states with different assignments to 
the leading/non-leading groups. 

In Fig. [5]Ja) we present the dependence of a on for 
low and high events for production by E~ beam 
events, to look for possible intrinsic charm effects p7] . 
Measurements of a are also available for D-meson pro- 
duction from the tt~ N experiments WA82 [11] (Si, Cu, 
W targets), E769 [H] (Be, Cu, Al, W), and WA92 [H] 
(Cu, W). We compare our £>- meson results from the tt" 
data to those from the other experiments in Fig. [51(b). 

As a systematic check, we performed the identical anal- 
ysis with and found good agreement with previous 
measurements [711^ for both proton and beams. De- 
tails will be presented in a forthcoming publication [29]. 
We looked for variations of a with any of the event selec- 
tion cuts. All changes were small compared with the sta- 
tistical uncertainty, indicating negligible systematic error 
from the cut selections. We also studied binning effects 
and found only small shifts, compatible with statistical 
uncertainties only. 



Table 4. a values for various grouping of charm species and/or 
beams. 

Average over xf and a 



all species 
all species 
all species 
all beams 
all beams 
all beams 

same as 
all beams 

same as 
At for E- 
A+ for E- 
D mesons 
D mesons 



and all beams 
for meson beams 
for baryon beams 
for charm particles 
for anti-charm particles 
for leading particles 
above, but excluding D^^s 
for non-leading particles 
above, but excluding D'^'s 
beam, p? < l.OGeVVc^ 
beam, p? > 1.3 GeVVc^ 
for ir" beam, pf < l.OGeVVc^ 
for TT- beam, p? > l.SGeVVc^ 



0.778 
0.850 
0.755 
0.763 
0.791 
0.814 
0.821 
0.747 
0.760 
0.894 
0.841 
0.836 
0.796 



±0.014 
± 0.028 
±0.016 
±0.021 
±0.019 
±0.021 
± 0.023 
±0.019 
± 0.034 
± 0.075 
±0.091 
± 0.045 
± 0.057 



5 Discussion and Conclusions 

As seen from the figures, all the measured values are com- 
patible with being independent of xp. Averaging over all 
our data, we obtain a — 0.778 ± 0.014, which is incom- 
patible with both usually suggested values of 2/3 and 1. 
Additional averages for different groupings are show in Ta- 
ble [4j Averaging separately over charm and anti-charm fi- 
nal states the a values show no difference. Separating into 
leading and non- leading production (Fig.|4]Jb)) we obtain 
a 2.3 cr difference, which vanishes due to the increased un- 
certainties when we exclude the Z^^'s in this comparison 
(Fig. HKc)). Separating the data into production by me- 
son beams and that by baryon beams (Fig. [4]Ja)), there 
is a difference in the a value averaged over all charm and 
anti-charm modes, corresponding to a 3 cr effect. 

For production of yl+ particles with a T,~ beam the 
behavior shown in Fig. OJa) seems to suggest a decrease 
for high Xp and , compared to D-mesons results shown 
in Fig. [5]^b). We note that _D-meson data has only a small 
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Fig. 4. Average a as a function oi xp for production by baryon {S~ , p) and meson (tt*) beams (a) and leading and non-leading 
(b, c) particles. The data points are slightly offset to avoid overlapping of the error bars. Reference a values of 2/3 and 1 are 
shown as dotted lines. The points at > 1 show the average assuming that a does not depend on xp- 
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Fig. 5. a for the production of with a E~ (a) and for D*, D*^ mesons with tv~ (b) beam, as function of xp for low 
(Pj < 1.0 GeV^/c^) and high (p^ > 1.3 GeV^/c^) transverse momentum. Our data points are slightly offset to avoid overlapping 
of the error bars. Reference a values of 2/3 and 1 are shown as dotted lines. Also shown (open symbols) are results |11II12|[T5] 
from other experiments, without separation in p^. 



contribution from events with pf > 1.3GeV^/c^, so no 
firm conclusion can be drawn. All distributions are con- 
sistent with no dependence on a;^- 

In summary, within our statistics there is no depen- 
dence of a on xf for any charm species for the inter- 
val 0.1 < Xf < 1-0. The average value of a for charm 
production by pion beams is Ofmeson = 0.850 ± 0.028. 
This is somewhat larger than the corresponding average 
c^baryon = 0.755 ± 0.016 for chami production by baryon 
beams {S~ , p). 
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